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Stress-Strain Relationships in Yarns Subjected to 
Rapid Impact Loading: 3, Effect of Wave 
Propagation ' 

Jack C. Smith, Frank L. McCrackin, and Herbert F. Schiefer 

The tensile behavior of a Hookean material, elongated by rapid impact at one end 
has been calculated, using a theory in which wave propagation is considered. As a result 
of these calculations, limits have been established on the applicability of a simpler theory, 
discussed in a preceding publication, in which wave propagation was neglected. 

1. Introduction 

When a filament of uniform cross section is clamped at each end and rapidly elongated, 
data for a stress-strain curve can be obtained by measuring the force at one end and the elonga- 
tion of the whole filament, both as functions of the time. In this case the assumption is made 
that stresses and strains in the filament are uniform along the length and can be represented 
by the stresses measured at an end and by the average elongations. This assumption is no 
longer valid when the rate at which the filament is lengthened approaches in order of magnitude 
the velocity of propagation of a tension wave along the filament. 

In previous work 2 this effect of wave propagation was disregarded. It is the purpose 
here to determine how previously obtained results must be reinterpreted as a result of the 
wave-propagation effect. 

2. Wave Theory 

Consider a filament lying along the negative x-axis, with one end at the origin and having 
the other end fixed at the point x=—L. At time £ = and subsequently, the end at x = is 
constrained to move with velocity v in the positive direction. Assume that the tensile behavior 
is in accordance with Hooke's law. Then, after impact, a strain pulse of magnitude e Q =Vo/c 
is propagated down the filament with a velocity c. When the pulse is reflected at the fixed 
end, the strain is increased to 2e . A similar increase takes place at each reflection. 

This behavior is shown in figure 1, in which the strains at several points along the filament 
are plotted as functions of the time. The strains are thus seen to vary discontinuously with 
time. These discontinuous increases in the strain become more pronounced as v is increased. 
If we require that e be less than 1 percent, in order that wave-propagation effects in most 
textile fibers may be safely disregarded, we must limit v to values less than 0.01 c, which in the 
case of high-tenacity nylon is approximately 25 m/sec. 

In the present experimental method a filament is terminated at x = —L by a mass nw, 
where n is a number and w is the mass of the filament. At #=0 the filament is terminated 
by a mass so large that its momentum will be practically unchanged by reaction forces from the 
filament and attached mass. At time £ — the head mass is given a velocity ?v The behavior 
of the system is recorded photographically. 3 The data of interest are the elongation of the 
filament and the reaction force applied to the filament at the tail mass, both as functions of 
time. From these data the experimental stress-strain curve is constructed. 



1 This work was sponsored by the Office of the Quartermaster General, Department of the Army, and was presented in part at the November 
1954 meeting of the Society of Rheology. 

2 Frank L. McCrackin, Herbert F. Schiefer, Jack C. Smith, and Walter K. Stone, Stress-strain relationships in yarns subjected to rapid 
impact loading: 2. Breaking velocities, strain energies, and theory neglecting wave propagation, J. Research NBS 54, 277 (1955) RP2590. 

3 Walter K. Stone, Herbert F. Schiefer, and George Fox, Stress-strain relationships in yarns subjected to rapid impact loading: 1. Equipment, 
testing procedure, and typical results, J. Research NBS 54, 269 (1955) RP2589. 
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Figure 1. Local strain versus time for a filament 
clamped at x=—L and rapidly elongated at x = 
at velocity v . 
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The theoretical behavior of this system, when Hooke's law applies, 4 is most easily calculated 
by finding suitable solutions of the wave equation 



dt 2 = 



dx 2 ' 



(1) 



In this equation, t is the time variable, and x denotes the position of a cross section of the 
filament relative to fixed coordinates, when the filament is in the unstrained state; u{x,t) 
denotes the distance the cross section has moved from its original position at x, and c is the 
velocity of propagation of an elastic-tension wave along the filament. This velocity is given 
by c—-y/E/p, where E is Young's modulus of elasticity for the unstrained filament, and p is 
the density of the unstrained filament material. 

These solutions must satisfy the boundary condition at the head : 



u(0, t) = v t. 
A suitable solution (for the unreflected wave) is given by 

u(x,t) = (t+~<0\ 

=*('+!) H>°) 

H>°> 



(2) 






_V0 

c 



(3) 



where u represents a wave of velocity c traveling along the filament in the minus z-direction, 
and e is the local strain, or increase in length per unit length, set up in the filament in the wake 
of the wave. It should be noted that both u and e are equal to zero for values of t<^\x\/c, i. e., 
no strain or displacement is present at a point in advance of the wave. 

When the wave arrives at the end of the filament it is reflected and travels back in the 
positive direction. The strain between the reflected wave front and the tail mass is represented 

* The stress-strain behavior of most textile fibers may be described by Hooke's law only for small strains. At larger strains, the solutions 
derived here provide only rough approximations to the behavior of an actual textile fiber but are adequate for the purposes of this theory. 

In order to take into account the significant reduction in cross-sectional area of the filament for the finite strains considered here, it is assumed 
that Hooke's law is obeyed in such a way that the product of Young's modulus and cross-sectional area remains constant. The small effect of mass 
motion associated with the lateral contraction has been neglected. 
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Figure 2. Strain-time curves for a filament loaded at 
the tail with a mass Ww and rapidly elongated at 
the head at velocity vq. 
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Figure 3. Stress-strain curves for a filament loaded 
with tail mass Ww and rapidly elongated at the head 
at velocity vq. 



as the sum of the strains caused by the incident and reflected waves. If the tail end of the 
filament were fixed, this strain would be 2e . However, as the tail mass can move, this strain, 
initially 2e Q at reflection, will dec revise with time. The rate of decrease depends upon the tail 
mass and upon the order of the reflection; i. e., whether it is the first or a subsequent reflection. 
Reflections also occur at the head mass. A different solution of the wave equation is required 
to describe the state of the filament after each reflection. The first 10 of these solutions are 
tabulated in appendix 1 . 

By using these solutions it is possible to calculate curves of average strain versus time, and 
local strain versus time. These curves are given in figure 2 for the case where the tail has a 
mass of lOw (n=\0). The curve of local strain, e(— Z, t), versus time is jagged, the decrease 
in e between successive reflections becoming greater with each reflection. 

As the force at the tail, which is found experimentally from the acceleration of the tail mass, 
is proportional to the strain at the tail, it is possible to construct a theoretical stress-strain 
curve of local tail stress (strain) versus average strain. Such a curve is given in figure 3. 
Experimental curves are similar to this curve in that they show dips in the vicinity of maximum 
stress. There is also a delay in the initial rise of stress by the time it takes for the tension wave 
to be propagated to the tail mass. 

3. Theory, Neglecting Wave Propagation 

A solution, applicable at low testing speeds, can be obtained by disregarding wave propaga- 
tion along the filament, and assuming the mass of the filament to be concentrated partly at the 
tail and partly at the head. The resulting solution, derived in appendix 2, is 



u= 



e L 0— V^' sin ~ 1 = 

L V^'J 



— =v« sin 



(4) 



(5) 



where d=(c/L)t, e =v /c, n , =n-\-f ) and / is the fraction of the filament mass assumed concen- 
trated at the tail. 

It is helpful to think of 6 as time expressed in L/c units. Likewise, one can think of [e/e ] as 
strain expressed in e units. Here e and c have only formal significance in order to put the 
solution in a form analogous to that obtained when wave propagation is considered. 
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4. Comparison of Wave Theory With Theory Neglecting Wave Propagation 

Plots of [e/eo] versus 6 for values of n' = 10, 10%, and 11 are shown in figure 4. Values of 
i/e have also been computed from the wave theory for the case 71=10, using the formula 



e_v4-u{-L,t) =d c 



60 



eoJ- 



(6) 



These points are plotted as circles and are seen to fit best to the curve drawn for n'=10%. 
This indicates that for a mass ratio of 10, the nonwave theory results agree most closely with 
those for the wave theory when half the mass of the filament is assumed concentrated at each 
end (/=*). 

According to the nonwave theory, the maximum average strain that can be obtained in a 
filament is given by 



L e 0j max 



(7) 



and this strain is obtained when m = (c/L)t m = (w/2)^/n+ \ or 



lm ^\ 7^ + 



(8) 



This suggests that the maximum average strain e/e as computed by the wave theory is approxi- 
mately equal to V^+i an d occurs for a value of 6 m approximately equal to t/2 V^+T- That 
this is a good approximation is demonstrated in table 1 . 
According to the nonwave theory (appendix 2) 



log V n = — \ log (n+i)+ log J- I " ade, 



(9) 



where v n is the lowest impact velocity to cause rupture in a filament to which a tail mass, nw, 
is attached, a is the tensile stress expressed as force divided by the cross-sectional area of the 
unstrained filament, and e T is the strain at which rupture occurs. 

Table 1. Values of maximum average strain and position of 
maximum average strain as computed by the wave theory 
and by the nonwave theory 



n 


Maximum average strain 


Position 


Wave theory 


Nonwave theory 


Wave theory 


Nonwave theory 




e/«o 


-yln+\ 


d m 


7r/2Vw+£ 


2 


1.61 


1.58 


2.39 


2.48 


5 


2.39 


2.34 


3.59 


3.68 


10 


3.26 


3.24 


5.10 


5.08 


25 


5.13 


5.05 


7.90 


7.93 



If the values of e r and a are independent of the rate of elongation, then log -%/(2/p) J crde 



will be a constant, and the plot of log v n versus log (n+|) will be a straight line of slope 



I 5 
2- 



5 At high rates of straining, the value of e T for an actual textile fiber is less than that obtained at conventional testing rates. Also the slope of 
the stress-strain curve is greater at higher rates of straining. However, for rates of straining of the same order of magnitude, these changes are 
small, and tend to cancel each other out. Consequently, the area under the stress-strain curve, or rupture energy density in the filament, should 
be essentially constant for rates of straining corresponding to impact velocities near the limiting breaking velocity. 
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When log {n-\-\) is zero, we have the extrapolated velocity 



When Hooke's law applies, we have 






(10) 

(11) 



According to the wave theory, -y/E/p is the velocity of propagation, c, of an elastic tension wave 
along the filament, and thus 

v b =ce r , 

the velocity just sufficient to cause immediate rupture at x=0. 

We call the velocity defined b^y (10) the 'limiting breaking velocity". The concept of a 
limiting breaking velocity has been discussed, in a previous paper (see footnote 2). 

As an example, consider a Hookean material for which € r =0.12, and c=2,500 m/sec. 



From (7), for the case of rupture, we have e T = € ->Jn-\- i = (v n /c) V^+ i or v n z 



Ce r 



V^+I 



(13) 



In figure 5 are plotted values of v n computed from these assumed values of c r and c. The 
result is a straight line giving an extrapolated value of 300 m/sec as the limiting breaking 
velocity. 

The breaking velocities have also been calculated on the basis of the more exact theory in 
which wave motion is considered and Hooke's law is assumed. Average strains are calculated 
from (6), and the breaking velocities are then computed from the formula 



C€ r 



L € J mux 



(14) 



These results are also plotted in figure 5. It is seen that these breaking velocities, for a 
Hookean material of constant e r , are lower (for small n's) than those predicted by the nonwave 
theory. 

In actuality the strain at a given point along the filament does not change continuously 
with time but increases suddenly by an amount e when one of the strain pulses passes through 
the point. Rupture then occurs whenever the local strain exceeds the rupture strain. 
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Figure 4. Strain-time curve for filament loaded with 
tail mass Ww, according to the wave theory, and 
according to the theory neglecting wave propagation. 




Figure 5. Breaking velocities for a hypothetical spec- 
imen, as computed by the wave theory, and by the 
nonwave theory. 
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Figure 6. — Local strain-time curves for a filament 
loaded with a tail-mass Ww. 



In figure 6 are plotted local strains at the head, tail, and midpoint for a filament loaded 
with a tail mass lOw (n=lO). Note here that maximum strain occurs at the head on the third 
reflection (0=6). This corresponds to a lowest breaking velocity of 68.6 m/sec. 

Increases in the local strain always occur first upon reflection at either the head or the tail. 
Increases in the local strain at an intervening point occur later, but these strains never exceed 
the head or tail strains in magnitude. For this reason, it is to be expected that a uniformly 
strong filament will break only at the ends. 

If the maximum values of local strains [e/e ]m&x are substituted into eq 14, the lowest 
breaking velocities, v n , for head and tail breaks can be computed. These breaking velocities, 
also plotted in figure 5, are seen to be the lowest of all. 

The curves of local strains versus time for an actual material, having approximately the 
properties chosen for this example, would not be as jagged as those shown in figure 6 for an 
ideal Hookean material. Thus, the velocities found experimentally for this actual material 
could be expected to lie between the bounds established in figure 5 by the curve of the velocities 
computed from average strains in the filament, and the curve of velocities computed from the 
extreme local strains at the head or tail of the filament. The breaking velocities computed by 
any of the methods considered are the same when n is greater than 50. It is thus possible, in 
this example, to determine a limiting breaking velocity by extrapolation of experimental data 
taken with n greater than 50, provided e r remains constant and Hooke's law applies. This 
criterion also applies, provided the energy density at rupture for the material remains nearly 
constant at the impact velocities used for determination of the limiting breaking velocity. 

In the example just considered, the breaking velocity for 71 = 50 is 42 m/sec. A tensile 
impact at this velocity produces a strain pulse of e =1.6 percent. This suggests a criterion 
applicable to any textile fiber, for which the energy density at rupture remains nearly constant 
at high strain rates. For such a fiber a limiting breaking velocity may be found by extrapola- 
tion of experimental data taken such that all v n used are less than 0.02 c, where c is the velocity 
of propagation of the tensile strain pulse along the sample for the material tested. 

5. Summary and Conclusions 

When the stress-strain characteristics of a textile fiber are measured at longitudinal impact 
rates exceeding 10 m/sec the effects of tension-wave propagation along the fiber must be con- 
sidered. These initial and reflected waves cause the local strains to vary with time in a step- 
wise fashion. The average strain in the filament, which is the quantity measured experi- 
mentally, does not adequately represent the state of strain under these conditions. As the 
measured stress is roughly proportional to the local strain at an end, experimental stress-time 
curves also vary in a stepwise fashion. Stress-strain curves likewise exhibit fluctuations. 

The tensile behavior of a Hookean material having a breaking elongation independent 
of testing speed has been calculated, using a theory in which wave propagation is considered. 
The results of these calculations have been used to modify the predictions of a simpler theory 
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that neglects wave propagation. Thus, it has been shown that a filament weighted at the tail 
and pulled rapidly at the head at constant speed will attain a maximum average elongation 
proportional to V^+l after the tension wave has been reflected (t/2)^ti-\- § times, where n is 
the ratio of the tail mass to the mass of the specimen. 

If logarithms of the lowest breaking velocities for each n are plotted versus log(?i+^), the 
theory neglecting wave propagation predicts that a straight line of slope — \ is obtained, 
provided that the rupture-energy density is nearly constant at high testing speeds. An extra- 
polated velocity (for log(?i+|) equals zero) is thus obtained, which is termed the "limiting 
breaking velocity'' of the material; i. e., the velocity above which the sample will always 
rupture at the head immediately upon impact. Both theories predict the same limiting break- 
ing velocity, provided extrapolation is made from data taken at impact speeds less than 2 
percent of the velocity of propagation of a tension-wave pulse along the material. If this limita- 
tion is observed, the theory neglecting wave propagation is valid, and extrapolation to the 
limiting breaking velocity may be made. 

6. Appendix 1. Wave-Theory Solutions 

In table 2, u gives the displacement due to the original incident wave, proceeding toward 
the tail mass. After reflection at the tail, the displacement, u, is computed from the sum of 
u , and U], the displacement of the reflected wave. Similarly, after the first head reflection, 
displacement is computed from the sum of u , u u and u 2 . Displacement due to each reflected 
wave is a function of a z parameter, which in turn is a function of x and t. Only positive values 
of these z parameters are allowed. Negative values represent position-time events in advance 
of the wave, where u=0. 

Velocity functions obtained by differentiating the displacement functions with respect to 
time are tabulated in table 3. Local-strain functions, obtained by differentiating the displace- 
ment functions with respect to the position coordinate along the filament, are given in table 4. 

The method of obtaining these solutions is demonstrated in the following derivation for u x . 

The initial conditions that must be satisfied for the first reflection at the tail are at 

x=—L, t=—) 2o=0, Zi = 0. 

Displacement at the tail just before reflection equals the displacement at the tail just after 
reflection, or u (0)=^o(0) + ^i(0). Hence U!(0) = 0. Velocity at the tail just before reflection 
equals the velocity at the tail just after reflection, or = # (0) + £>i(0). Hence v l (Q)=—v . In 
addition, we have the boundary condition at x=—L, 



'5i 2rw ' *dx 

or 

w u l {z)=EA- 



nw^^[uo(z Q )+u 1 (zi)] x =. L =EA^ z [uo(zo) J rU l (z l )] x ^. L 



EA^- x uM)\ x 



or 

d 2 _ A v EA d 

dz z c c dz 



where z=t — L/c, and A is the cross-sectional area of the unstrained filament. If we let 

nwc npLc L 

we obtain 

cPui . 1 dui_v 
dz 2 r dz r 
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The solution of this equation that satisfies the initial conditions is 

If we replace (t—L/c) by z 1 = t—2(L/c)—x/c, we obtain the required solution for %. 



Table 2. Displacement functions 

u = u + Ui + u 2 + Uz J r .... 



U = VoZq 

Ui=v Zi— 2v t\ l—e~ T 
u 2 =—VoZ 2 +2vot\ l—e~ T 



y where r=n : 



u z =— v z z +4vqt\ 1— (l+^Je * 

u 4 =v z A — 4v t\ 1— (l+— )e~ T 

u 5 =VoZ 5 -6voT\l-U+-^+--^)e *J 

u 6 =-v z 6 +6vot[i-(i+ 2 3 7+| 7)^] 

u 7 ^-V z 7 +Svor^l-(^+^+l^y-^ 

u 8 =VoZ 8 -Sv T^l-(^+^+^y-^ 

*« T 1 A ,409, 4^1 4 2j . 2 gj\ -!?1 

.in l~i A , 4 210,4 2? 4 2* 2 4A _£f| 
Wlo== -_^ 10+ 10, or ^l-^l+- _+- _-_ _+_ v )e+j 



Table 3. Velocity functions 



_ du_ duo . dwi , dw 2 1 



v =v (for 2 >0) 

z?2=— #o+2#oe T 



^3=— ^o+4^o— e' 

r 



5.-5 



«b=*+| 



c c 

. 2L.x 

z 2 =t — 

c c 

4i x 

z z =t 

c c 

. 4L x 

z±=t \— 

c c 

Zb ~ c c 
. 6i . x 

Zl ~ c c 
. 8Z . x 



Z9 = t- 



10L g 



. 10Z . x 

2io=* h- 

c c 
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Table 3. Velocity functions — Continued 



v 4 =v — 4?; — e T 

T 






. a (\ 2 2 6| 2 gj\ _£5 

in /l 4 2 9 , 8 2? 4 z| , 2 4\ Ji 
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Table 4. Local strain functions 

du._duo.du! _ 

c — >r~ — -5: — r ^ — r • • • — €oi-€ii- . . . 
ox ox dx 



t0 =|» (for 2o >0) 



_22 

C2= — eo+2e e r 



€ 3 =e — 4e — e~ 

T 



€ 4 =€ — 4e — 6 T 

T 

e5= _ 6o+6eo ^---_ + -^e 

. „ /l 2 2 6 , 2 2 2 A _£• 

«6 = — «0+O€ l g — g — +g ^2 K T 

€ 7 =«o-8«o^--^+3^jc r 
/z 8 2§ , 1 2Jj\ _a 

/l 4 2 9 8 2 2 9 4 2|,2 2 9 \ * 

e9 __ €o+ l 0eo ^-_--- 1 -- ? _- ?+ --je 

, - *+in, I" 1 42 >»i 82 ?o 4 2? 2 2t \ » 
*•" ~ €0+ ! ° £o V5~5 T+5 "5 ^ + 15 V ) 6 
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Appendix 2. Solution, Neglecting Wave Propagation 

Consider a filament of length L, cross section A, and density p and having a mass n' pAL 
attached at the tail, where n r is equal to n-\-j, n is the ratio of the tail mass to the mass of the 
filament, and / is the fraction of the mass of the filament assumed concentrated at the tail. It 
is shown in section 4 that the best value for /is |. 

At zero time and subsequently, let the head end of the filament be constrained to move with 
velocity v . The equation of motion is 

n' P AL d ^=F, (15) 

where u is the displacement of the tail mass, and F is the accelerating force transmitted through 
the filament to the tail mass. 

The strain in the filament, e, is given by 

vpt—u 
e= -£~- (16) 

When Hooke's law applies, 

F=AEe, (17) 

where E is Young's modulus for the filament material. Substituting (16) and (17) into (15), 
we obtain 

d 2 u 

-jp+uioU= <$><$, (18) 



where 



w „=4 -L (19) 



and 

IE 



V 



(20) 
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A solution of (18) satisfying the boundary conditions, u=0 and du/dt=0 at £— 0, is 

u=v t sin coot • (21) 

If we let t=(L/c)0, and e =Vo/c, we obtain 

u = e L 6—4n' sin — (22) 

a 

e=€oV^' sin -=• (23) 

In a previous paper (see footnote 2) the following equation was derived: 



where v n is the lowest impact velocity sufficient to cause rupture in a filament to which a tail 
mass n times that of the filament and a head mass m times that of the filament, are attached. 
For comparison with the work discussed here, we let ra— > «> , obtaining 



log v n =-\ log (n+h)+\og ^ J' r ade. (25) 

Washington, December 31, 1954. 
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